0^ 



O 



Binding Energies and Scattering Observables in the ^Hes Atomic System* 

A. K. Motovilovt, W. Sandhas 

Physikalisches Institut der Universitdt Bonn, Endenicher Allee 11-13, D-53115 Bonn, Germany 

S. A. Sofianos 

Physics Department, University of South Africa, P.O.Box 392, Pretoria 0003, South Africa 



E. A. Kolganova 

Joint Institute for Nuclear Research, Dubna, 141980, Russia 
(November 8, 1999) 
^ . 
O ■ The *He3 bound states and the scattering of a *He atom off a ''He dimer at ultra-low energies 

1/'^ ' are investigated using a hard-core version of the Faddeev differential equations. Various realistic 

He- He interactions were employed, amomg them the LM2M2 potential by Aziz and Slaman and 
the recent TTY potential by Tang, Toennies and Yiu. The ground state and the excited (Efimov) 
state obtained are compared with other results. The scattering lengths and the atom-diatom phase 
shifts were calculated for center of mass energies up to 2.45 mK. It was found that the LM2M2 and 
r^ ' TTY potentials, although of quite different structure, give practically the same bound-state and 

^1^' scattering results. 

PACS numbers: 02.60.Nm, 21.45. -fv, 34.40.-m, 36.40. -Hd 



I. INTRODUCTION 



CJ 

^ . Small ^He clusters (in particular dimers and trimers) are of fundamental interest in various fields 

^ ' of physical chemistry and molecular physics. Studies of these clusters represent an important step 

Ph towards understanding the properties of helium liquid drops, super-fluidity in ^He films, the Bose- 

Einstein condensation etc. (see, for instance, Refs. |]T]-0)- Besides, the helium trimer is probably a 

CN ,: unique molecular system where a direct manifestation of the Efimov effect |^ can be observed since 

vQ ] the binding energy e^ of the '^He dimer is extremely small. 

T-H I The ''He trimer belongs to the three-body systems whose theoretical treatment is quite difficult, 

^ ■ first, due to its Efimov nature and, second, because of the hard-core properties of the inter-atomic 

,__! . He -He interaction ^-^. At the same time the problem of three helium atoms can be considered as 

0^ [ an example of an ideal three-body quantum problem since ^He atoms are identical neutral bosons 

^ ' and, thus, their handling is not complicated by spin, isospin, or Coulomb considerations. 

Q There is a great number of experimental and theoretical studies of ''He clusters. However, most 

^ of the theoretical investigations consist merely in computing the ground-state energy and are based 



>>^, on variational methods ||TO|-[T^, on Hyperspherical Harmonics expansion methods in configuration 
Oh. 



space ||16|JT7|, and on integral equations in momentum space [18,19|. We further note that the 



results of Ref. 20 were based on a direct solution of the two-dimensional Faddeev differential 
• 1— I . equations in configuration space, while recently binding-energy results were obtained by using the 
r> ■ three-dimensional Faddeev differential equations in the total-angular-momentum representation [^ . 



c^ A qualitative treatment of the ^Hea system in framework of an effective field theory is presented in 
Ref. [0. 
In Refs. p3|,p!6|,p!9|,p3| it was pointed out that the excited state of the ^He trimer is an Efimov state 



||. In these works the HFDHE2 0, HFD-B §, and LM2M2 g] versions of the ^He-^He potentials 
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by Aziz and co-workers were employed. The essential property of this state is that it disappears 
when the inter-atomic potential is increased by a factor A ~ 1.2. And vice versa, when A slightly 
decreases (no more than 2 %), a second excited state appears in the trimer |]T6|,p!9|. It is just such a 



non-standard behavior of the excited-state energies which points at their Efimov nature. Another 
proof of such a nature of the excited state of the trimer, based on a scaling consideration, is discussed 
in 1^1 . The resonance mechanism of formation and disappearance of the Efimov levels in the ^Hes 
system has been studied in Refs. |25|j2r 



The general atom-diatom collision problem has been addressed by various authors, and we refer 
the interested reader to the review articles [^] and [^ . The collision dynamics at thermal energies 
of the H-I-H2 system and the existence of resonances were discussed in |^ using the Faddeev integral 
equations in momentum space. Finally, the problem of existence of ^He n-mers and their relation to 
the Bose-Einstein condensation in Hell was discussed in Refs. [3D, JT|. From the experimental studies 
we mention those of Refs. [^^ -|37|| where molecular clusters, consisting of a small number of noble 
gas atoms, were investigated. 

In contrast to the bulk of theoretical investigations devoted to the binding energies of the ^He 
trimer, scattering processes found comparatively little attention. In Ref. [|18| the characteristics of 
the He-He2 scattering at zero energy were studied, while the recombination rate of the reaction 
(1 + 1 + 1 -^ 2 + 1) was estimated in 13^. Recently, the phase shifts of the He-He2 elastic scattering 



and breakup amplitudes at ultra-low energies have also been calculated p3| , p9| , ^ |. 

The difficulty in computing excited states and scattering observables in the ^Hes system is mainly 
due to two reasons. First, the low energy e^ of the dimer makes it necessary to consider very 
large domains in configuration space with a characteristic size of hundreds of Angstroems. Second, 
the strong repulsion of the He-He interaction at short distances produces large numerical errors. 
In the present paper, which is an extension of our studies for ^Hea p3|,p5|,p6|,|39|,^ , we employed 



the mathematically rigorous three-body Boundary Condition Model (BCM) of Refs. [p!| , ^ to the 
above-mentioned problems. 

As compared to |2^J2^J2^J3^,^ we employ, in the present work, the refined He-He interatomic 
potentials LM2M2 by Aziz and Slaman |^ , and TTY by Tang, Toennies, and Yiu [^ . Our numerical 
methods have been substantially improved, and this allowed us to use considerably larger grids 
achieving, thus, a better accuracy. Furthermore, due to much better computing facilities more 
partial waves could be taken into account. 

This paper is organized as follows. In Sec. || we review the three-body bound and scattering state 
formalism for hard-core interactions. In Sec. |T| we describe its application to the system of three 
^He atoms and present our numerical results. Our conclusions are summarized in Sec. |^. Finally 
in the Appendix we give details of the potentials used. 



II. FORMALISM 



A detailed analysis of the general boundary- value problem, the derivation of the asymptotic bound- 
ary conditions for scattering states and other boundary-value formulations, can be found in Refs. 
43| , |4^ . In this work we employ a hard-core version of the BCM P2| , |i5[] developed in [p] , ^ (for 
details see Ref. [pH]). Therefore, in what follows we shall only outline the formalism and present its 



main characteristics. 

In describing the three-body system we use the standard Jacobi coordinates XQ,,yQ,, a 
expressed in terms of the position vectors of the particles r^ and their masses m^. 



1,2,3, 



nl/2 




m/3 + m^ 

where (a, /?, 7) stands for a cyclic permutation of the indices (1, 2, 3). 

In the so-called hard-core potential model one requires that the wave function vanishes when 
the particles approach each other at a certain distance r = c. This requirement is equivalent to 
the introduction of an infinitely strong repulsion between the particles at distances r < c. Such a 
replacement of the repulsive short-range part of the potential by a hard-core interaction turns out 
to be a very efficient way to suppress inaccuracies at short distances. One can then show that the 
Faddeev components satisfy the following system of differential equations 

/3^a (1) 

(-Ax-^)$a(X) = 0, |X,| <Ca. 



where X 




1,2,3 and Cn is the hard-core radius in the channel a. 



Outside the core the components $a still provide the total wave function \E', 



while in the interior region we have 



(3=1 



(3=1 



In practice, one can replace the latter strong condition by a more weak one |41,46 



Y.MX) 



(3=1 



a 



1,2,3 



(2) 



|Xc«| = 



which requires the sum of $a(X) to be zero only at the radius Cq. 

The numerical advantage of our approach is already obvious from the structure of Eqs. (^. When 
a potential with a strong repulsive core is replaced by the hard-core model, one approximates inside 
the core domains only the Laplacian Ax instead of the sum of the Laplacian and the huge repulsive 
term. In this way a much better numerical approximation can be achieved. 

In the present investigation we apply the formalism to the ^He three-atomic system with total 
angular momentum L = 0. The partial- wave version of the equations (|^) for a system of three 
identical bosons with L = reads P7| , |i5[ | 



d^ 



d^ 



dx^ dy" 



+ l{l + l) 



1 1 

y 



x^ 



-E 



'^i{x,y) 



-V{x)'^i{x,y), x> c 



0, 



X < c. 



(3) 



Here, x, y are the absolute values of the Jacobi variables and c is the core size which is the same for 
all three two-body subsystems. The angular momentum / corresponds to a dimer subsystem and a 
complementary atom. For a three-boson system in an S'-state / can only be even, / = 0, 2, 4, . . . . The 



potential V{x) is assumed to be central and the same for all partial waves /. The function \E';(x,?/) 
is related to the partial- wave Faddeev components ^i{x,y) by 

'ifi{x,y) = <l>i{x,y) + J2 ^V hw {x,y,r]) $;/ (x', y') (4) 

where 



, ; 1 2 So V 3 / /So 1 2 V 3 



with ?7 = X ■ y . Expressions for the kernels /i;;/ can be found in p3| , ^^ . It should be noted that 
these kernels depend only on the hyperangles 

6 = arctan — and 6' = arctan — 
X x' 

and not on the hyperradius 

p = \/xP^+y^ = a/o;'^ + y'"^. 
The functions ^i{x,y) satisfy the boundary conditions 

The partial-wave version of the hard-core boundary condition (j^) reads 

r+l 

<^i{c,y) + J2 dr]hu,{c,y,7])Mx',y') = (6) 

I' -^"1 

requiring the wave function '^i[x,y) to be zero at the core boundary x = c. Furthermore, one can 
show that, in general, the condition (|]), like the condition (H), causes also the wave functions (|) to 
vanish inside the core domains. For the bound-state problem one requires that the functions $/(a;, y) 
are square integrable in the quadrant x > 0, y > 0. 
The asymptotic condition for the helium trimer scattering states reads 

^i{x,y) = 6ioipd{x) exp{iy/Et - eay) [ao + o (y"^/^)] 

^ exp(,^p) , (7) 

as p — i> oo and/or y ^ oo. Here we use the fact that the helium dimer bound state exists only for 
/ = 0. td stands for the dimer energy while il^d{x) denotes the dimer wave function which is assumed 
to be zero within the core, i. e., ipd{x) = for x < c. 

The coefficients ao and Ai{6) describe the contributions of the (2 + 1) and (1 + 1 + 1) channels to 
$/, respectively. Both the trimer binding energy Et and the difference Et — e^ in (|^ are negative 
which means that for any 6 the function ^i{x,y) decreases exponentially as p ^ oo . 

The asymptotic boundary condition for the partial-wave Faddeev components of the (2 + 1 — > 
2 + 1; 1 + 1 + 1) scattering wave function reads, as p ^ oo and/or y ^ oo, 

^i{x,y;p) = diQipdix) {sm{py) + exp{ipy) [aoip) + o {y'^^^)] } 



^ej^pl^VEpl (8) 

\fP 



where p is the relative momentum conjugate to the variable y, E is the scattering energy given by 
E = e^ + p^, and ao(p) is the elastic scattering amplitude. The functions Ai{E,6) provide us for 
E > {] with the corresponding partial-wave breakup amplitudes. 
The helium-atom helium-dimer scattering length l^c is given by 

2 p 

p^ ^ 

while the S'-state elastic scattering phase shifts 5q{p) are given by 

5o(p) = ^ImlnSo(p) (9) 

where So(p) = 1 + 2iao(p) is the (2 + 1^2 + 1) partial-wave component of the scattering matrix. 

III. RESULTS 

We employed the Faddeev equations (D and the hard-core boundary condition (^ to calculate 
the binding energies of the helium trimer and the ultra-low energy phase shifts of the helium atom 
scattered off the helium diatomic molecule. As He-He interaction we used three versions of the 
semi-empirical potentials of Aziz and collaborators, namely HFDHE2 P], HFD-B 0, and the newer 
version LM2M2 |^. Further, we employed the latest theoretical He-He potential TTY of Tang 
et al. 0. These potentials are given in the Appendix. In our calculations we used the value 
fi^/ra. = 12.12 KA^. All the potentials considered produce a weakly bound dimer state. The energy 
ed of this state together with the He-He atomic scattering length ilc are given in Table |. It is 
interesting to note that the latest potentials LM2M2 and TTY give practically the same scattering 
length 4c and dimer energy e^. 

A detailed description of our numerical method has been given in Ref. |23|. Therefore, we outline 



here only the main steps of the computational scheme employed to solve the boundary- value problems 
(|), (|^), (|) and (0) or (H). First, we note that the grid for the finite-difference approximation of 
the polar coordinates p and 6 is chosen such that the points of intersection of the arcs p = pi, 
i = 1,2, ... ,Np and the rays 9 = 9j, j = 1,2, ... ,Ng with the core boundary x = c constitute 
the knots. The value of c was fixed to be such that any further decrease of it did not appreciably 
influence the dimer binding energy e^ and the energy of the trimer ground state E^ . In our previous 
work [^,^,0 c was chosen as 0.7 A. In the present work, however, we choose c = 1.0 A. This value 
of c provides a dimer bound state e^ which is stable within six figures and a trimer ground-state 
energy E^ stable within three figures. The pi are chosen according to the formulas 



Pi 



M'^ + 1 



t = l,2,...,N^''\ 



P^+Ni^^ = V ^' + ^^'' ^ = 1, 2, . . . , AT, - N^P\ 
where Nc stands for the number of arcs inside the domain p < c and 



yi = f{n)\ pI^-c^, n 



' N. - n!^'^ ' 



The nonlinear monotonously increasing function /(t), < r < 1, satisfying the conditions /(O) =0 
and /(I) = 1, is chosen according to 

f(^\ ^ / "07- , T e [0, To] 

^^ ' [aiT + T'' , rG(ro,l] ' 

The values of ao, a^ > 0, and ai, ai > 0, are determined via tq and u from the continuity condition 
for /(r) and its derivative at the point tq. In the present investigation we took values of tq within 
0.15 and 0.2. The value of the power v depends on the cutoff radius pmax = pNp = 200-1000 A, its 
range being within 3.4 and 4 in the present calculations. 
The knots 9j at j = 1,2, . . . , Np — Nc are taken according to 6j = arctan(?/j/c) with the remaining 

knots 9j, j = Np — Nc + 1, . . . , A^g, being chosen equi distantly. Such a choice is required by the need 
of having a higher density of points in the domain where the functions ^i{x, y; z) change most rapidly, 
i. e., for small values of p and/or x. In this work, we used grids of the dimension Nq = Np = 500-800 
while the above number Nc and the number Ng — {Np — Nc ) of knots in 9 lying in the last arc 
inside the core domain was chosen equal to 2-5. 

Since we consider identical bosons only the components $; corresponding to even / differ from 
zero. Thus, the number of equations to be solved is A^e = ^max/2 + 1 where /max is the maximal even 
partial wave. The finite-difference approximation of the A'"e equations (|^) reduces the problem to a 
system of N^NoNp linear algebraic equations. The finite-difference equations corresponding to the 
arc i = Np include initially the values of the unknown functions ^i{x, y; z) from the arc i = Np + 1. 
To eliminate them, we express these values through the values of ^i{x, y; z) on the arcs i = Np and 
i = Np — 1 hj using the asymptotic formulas (^ or (|^) in the manner described in the final part 



of Appendix A of Ref. p3|. In ||2^, however, this approach was used for computing the binding 
energies only, while in the present work this method is extended also to the scattering problem. The 
matrix of the resulting system of equations has a block three-diagonal form. Every block has the 
dimension NqNq x NeNo and consists of the coefficients standing at unknown values of the Faddeev 
components in the grid knots belonging to a certain arc p = pi. The main diagonal of the matrix 
consists of Np such blocks. 

In this work we solve the block three-diagonal algebraic system on the basis of the matrix sweep 
method |^^. The use of this method makes it possible to avoid writing the matrix on the hard drive 



of the computer. Besides, the matrix sweep method reduces the computer time required by almost 



one order of magnitude as compared to p3|j39|j40[| . 



Our results for the trimer ground-state energy E^ , as well as the results obtained by other authors, 
are presented in Table |I[ It should be noted that most of the contribution to the ground-state energy 
stems from the / = and / = 2 partial-wave components, the latter being slightly more than 30 %, 
and is approximately the same for all potentials used. The contribution from the / = 4 partial wave 
is of the order of 3-4 % (cf. 0). 

It is well known that the excited state of the ^He trimer is an Efimov state |l6| , p!9| , |23| , p5| , |26| . The 
results obtained for its energy E^ , as well as the results found in the literature, are presented in 
Table |TT[ To illustrate the convergence of our results we show in Table |rV| the dependence of El 
on the grid parameters using the TTY potential. It is seen that the / = partial-wave component 
contributes about 71 % to the excited-state binding energy. The contribution to E^ from the / = 2 
component is about 25-26 % and from / = 4 within 3-4 %. These values are similar to the ones for 
the ground state. 

Besides the binding-energy calculations, we also performed calculations for a helium atom scattered 
off a helium dimer for L = 0. For this we used the asymptotic boundary conditions (H). The 



scattering lengths of the colhsion of the He atom on the He dimer obtained for the HFD-B, LM2M2 
and TTY potentials are presented in Table 0. As compared to p3| the present calcualtion is 



essentially improved (the result 4c = 145±5A for HFD-B with /max = 2 was obtained in [23| with 
a much smaller grid). Within the accuracy of our calculations, the scattering lengths provided by 
the LM2M2 and TTY potentials, like the energies of the excited state, are exactly the same. This is 
of no surprise since the two potentials produce practically the same two-body binding energies and 
scattering lengths. 

The phase-shift results obtained for the HFD-B, LM2M2 and TTY potentials are given in Tables 
0, [VT], and [VITI| . For the HFD-B and TTY potentials they are plotted in Fig. |I]. Note that for the 
phase shifts we use the normalization required by the Levinson theorem |^, Sl{0) — Sl{oo) = nn, 



where n is the number of the trimer bound states. 

Incident energies below and above the breakup threshold were considered, i.e., for the (2-1-1 — > 
2-1-1) and (2-1-1 — > 1 -|- 1 -|- 1) processes. It was found that after transformation to the laboratory 

system the phases 5^'"'' for the potentials HFD-B, LM2M2 and TTY for different values of /max are 
practically the same, especially those for LM2M2 and TTY. The difference between the phase shifts 
6q and 6^ is only about 0.5 %. 

It is interesting to compare the values obtained for the He-He2 scattering lengths 4c with the 
corresponding inverse wave numbers x~^ for the trimer excited-state energies. The values of x^^ 



where x = 2^J [td — El )/3, with both the Ei and e^ being given in A ^, are also presented in 
Table |V|. It is seen that these numbers are about 1.3-1.7 times smaller than the respective ^He-atom 
^He-dimer scattering lengths. The situation differs completely from the ^He two-atomic scattering 
problem, where the inverse wave numbers (x^^^)" = je^l"^/^ are in a rather good agreement with 
the ^He-^He scattering lengths (see Table |). Such significant differences between 4c and x~^ in the 
case of the '^He three-atomic system can be attributed to the Efimov nature of the excited state of 
the trimer, which implies that the effective range tq for the interaction between the ^He atom and 
the "^He dimer is very large as compared to the ^He diatomic problem. 

IV. CONCLUSIONS 

In this paper we employed a hard-core variant of the Boundary Condition Model which, unlike 
some competing methods, is exact and ideally suited for three-body calculations with two-body 
interactions characterized by a highly repulsive core. The method is relevant not only for bound 
states but also for scattering processes. There is, however, a price to be paid for the exact treatment 
of the system. Higher partial waves, beyond /max = 4, were hard to be incorporated within the 
computing facilities at our disposal. 

The results for the ground-state energy of the '^He trimer obtained for all four '^He-'^He potentials 
considered compare favorably with alternative results in the literature. Furthermore, the treatment 
of the excited state, interpreted as an Efimov state, clearly demonstrates the reliability of our method 
in three-body bound state calculations with hard-core potentials. In addition to the binding energy, 
the formalism has been successfully used to calculate scattering lengths and ultra-low-energy phase 
shifts of the ^He atom scattered off the ^He dimer. 

In general the hard-core inter-atomic potential together with other characteristics of the system, 
makes calculations extremely tedious and numerically unstable. This is not the case in our formalism 
where the hard core is taken into account from the very beginning in a mathematically rigorous way. 
Thus, the formalism paves the way to study various ultra-cold three-atom systems, and to calculate 
important quantities as cross-sections, recombination rates, etc. 
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APPENDIX: THE POTENTIALS USED 

The general structure of the realistic semi-empirical potentials HFDHE2 |^ and HFD-B 0] devel- 
oped by Aziz et al. is 



V{x)=eV,{0 
where C = a^/^m and the term Vfe(C) reads 



(Al) 



H(C)=Aexp(-«C + /9C' 






no 



X 



is expressed in the same length units as r^ (A in the present case). The function -F(C) is given by 



HO 



exp[-{D/C-l)f, ifC<^ 



1, iiC> D. 

In addition to the term Vb{0 ^^e LM2M2 potential contains a term V^(C), 

y(r)=e{H(C) + K(C)}, 
where 

"27r(C-Ci) vr 



(A2) 



VaiO 



Aa { sin 



C2-C1 
0, 



1 



Ci < C < C2 

C^[Ci,C2]. 



The parameters for the HFDHE2, HFD-B and LM2M2 potentials are given in Table ^. 
The form of the theoretical He-He potential TTY is taken from 0. This potential reads 

V{x) = A %^{x) + Vdisp(x)] 

where x stands for the distance between the '^He atoms given in atomic length units. (Following 
in converting the length units we used la.u.= 0.52917 A.) The function V^^ has the form 

V^^{x) = D xP exp{-2px) 



with p = — — 1 , while the function Vdisp reads 

2p 



N 



Vdisp (a;) 



7 , C2n f: 



2n[X X 



-2n 



n=3 



The coefficients C2„ are calculated via the recurrency relation 



<-'2n — -p. L^2n- 



a 



2n-4 



and the functions f2n are given by 



2n 



/2„(x) = 1 - exp(-6x) y^ 



A;=0 






where 



h-^ 



h{x) = 2j3- 
The parameters of the TTY potential are given in Table ^. 



1 

X 



[1 

[2; 
[3; 
[4; 
[s; 
[6; 
[r 

[9; 
[lo; 
[11 

[12 
[13 
[14 
[15 

[16 
[17 
[18 
[19 
[20 

[21 
[22' 
[23^ 
[24' 
[25^ 
[26^ 
[27 
[28' 
[29 
[30 

[31 
[32' 
[33 

[34; 
[35; 



McConville, J. Chem. Phys. 79, 4330 (1979). 
61, 1487 (1987). 



J. P. Toennies and K. Winkelmann, J. Chem. Phys. 66, 3965 (1977). 

M. V. Rama Krishna and K. B. Whaley, Phys. Rev. Lett. 64, 1126 (1990). 

K. K. Lehman and G. Scoles, Science 279, 2065 (1998). 

S. Grebenev, J. P. Toennies, and A. F. Vilesov, Science 279, 2083 (1998). 

V. Efimov, Nucl. Phys. A 210, 157 (1973). 

R. A. Aziz, V. P. S. Nain, J. S. Carley, W. L. Taylor, and G. T. 

R. A. Aziz, F. R. W. McCourt, and C. C. K. Wong, Mol. Phys. 

R. A. Aziz and M. J. Slaman, J. Chem. Phys. 94, 8047 (1991). 

K. T. Tang, J. P. Toennies, and C. L. Yiu, Phys. Rev. Lett. 74, 1546 (1995). 

S. W. Rick, D. L. Lynch, and J. D. Doll, J. Chem. Phys. 95, 3506 (1991). 

V. R. Pandharipande, J. G. Zabolitzky, S. C. Pieper, R. B. Wiringa, and U. Helmbrecht, Phys. Rev. Lett., 50, 1676 (1983). 

R. N. Barnett and K. B. Whaley, Phys. Rev. A 47, 4082 (1993). 

M. Lewerenz, J. Chem. Phys. 106, 4596 (1997). 

R. Guardiola, M. Portesi, and J. Navarro, LANL E-print 3hysics/9904037 . 

T. Gonzalez-Lezana, J. Rubayo-Soneira, S. Miret-Artes, F. A. Gianturco, G. Delgado-Barrio, and P. Villareal, Phys. Rev. 

Lett. 82, 1648 (1999). 

B. D. Esry, C. D. Lin, and C. H. Greene, Phys. Rev. A 54, 394 (1996). 

E. Nielsen, D. V. Fedorov, and A. S. Jensen, J. Phys. B 31, 4085 (1998) (LANL E-print |physics/98O6O20| ). 

S. Nakaichi-Maeda and T. K. Lim, Phys. Rev A 28, 692 (1983). 

Th. Cornelius and W. Glockle, J. Chem. Phys. 85, 3906 (1986). 

J. CarboneU, C. Gignoux, and S. P. Merkuri ev, Few-Body Sy stems 15, 15 (1993). 



physics/980802 



V. Roudnev and S. Yakovlev, LANL E-print ^hysics/991003C 

P. F. Bedaque, H.-W. Hammer, and U. van Kolck, Nucl. Phys. A 646, 444 (1999) (LANL E-print |nucl-th/9811046|) . 

E. A. Kolganova, A. K. Motovilov, and S. A. Sofianos, J. Phys. B. 31, 1279 (1998) (LANL E-print ^lysics/9612012| ) 

T. Frederico, L. Tomio, A. Delfino, and A. E. A. Amorim, Phys. Rev. A 60, R9 (1999). 

A. K. Motovilov and E. A. Kolganova, Few-Body Systems Suppl. 10, 75 (1999) (LANL E-print [3hvsics/981000i 

E. A. Kolganova and A. K. Motovilov, Phys. Atom. Nucl. 62 No. 7, 1179 (1999) (LANL E-print 
D. A. Micha, Nucl. Phys. A 353, 309 (1981). 
A. Kuppermann, Nucl. Phys. A 353, 287 (1981). 
Z. C. Kuruoglu and D. A. Micha, J. Chem. Phys. 80, 4262 (1980). 
H. B. Ghassib and G. V. Chester, J. Chem. Phys. 82, 585 (1984). 
N. H. March, J. Chem. Phys. 82, 587 (1984). 

F. Luo, G. C. McBane, G. Kim, C. F. Giese, and W. R. Gentry, J. Chem. Phys. 98, 3564 (1993) 
F. Luo, C. F. Giese, and W. R. Gentry, J. Chem. Phys. 104, 1151 (1996). 
W. Schollkopf and J. P. Toennies, Science 266, 1345 (1994). 
U. Buck and H. Meyer, J. Chem. Phys. 84, 4854 (1986). 



[36] O. Echt, K. Sattler, and E. Recknagel, Phys. Rev. Lett. 47, 1121 (1981). 

[37] W. Schollkopf and J. R Toennies, J. Chem. Phys. 104, 1155 (1996). 

[38] P. O. Fedichev, M. W. Reynolds, and G. V. Shlyapnikov, Phys. Rev. Lett. 77, 2921 (1996). 

[39] E. A. Kolganova, A. K. Motovilov and S. A. Sofianos, Phys. Rev. A. 56, R1686 (1997) (LANL E-print |physics/98020ie|] 



[40] A. K. Motovilov, S. A. Sofianos, and E. A. Kolganova, Chem. Phys. Lett. 275, 168 (1997) (LANL E-print |physics/9709037[ ) 

[41] S. P. Merkuriev and A. K. Motovilov, Lett. Math. Phys. 7, 497 (1983). 

[42] S. P. Merkuriev, A. K. Motovilov, and S. L. Yakovlev, Theor. Math. Phys. 94, 306 (1993). 

[43] A. K. Motovilov, Three-body quantum problem in the boundary-condition model (PhD thesis (in Russian), Leningrad State 

University, Leningrad, 1984). 
[44] A. A. Kvitsinsky, Yu. A. Kuperin, S. P. Merkuriev, A. K. Motovilov, and S. L. Yakovlev, Sov. J. Part. Nucl. 17, 113 

(1986). 
[45] V. N. Efimov and H. Schulz, Sov. J. Part. Nucl. 7, 349 (1976). 
[46] A. K. Motovilov, Vestnik Leningradskogo Universiteta, 22, 76 (1983). 
[47] L. D. Faddeev and S. P. Merkuriev, Quantum scattering theory for several particle systems (Doderecht: Kluwer Academic 

Publishers, 1993). 
[48] S. P. Merkuriev, C. Gignoux, and A. Laverne, Ann. Phys. (N.Y.) 99, 30 (1976). 
[49] A. A. Samarsky: Theory of difference schemes (in Russian) (Nauka, Moscow, 1977). 
[50] N. Levinson, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 25, 9 (1949). 



10 



TABLE I. 


Dimer energies e^, 


inverse wave lengths 1/x^^ 


, and "He-^He scattering lengths ^^c^ 


for the potentials used. 


Potential 




Ei (mK) 


i/x(2) (A) ej (A) 


Potential 


Ei (mK) 


i/x(2) (A) ej (A) 


HFDHE2 
HFD-B 




-0.83012 
-1.68541 


120.83 124.65 
84.80 88.50 


LM2M2 
TTY 


1.30348 
1.30962 


96.43 100.23 
96.20 100.01 



TABLE II. Ground state energy E\ results for the helium trimer. The (absolute) values of El are given in K. The grid 



(0) 



parameters used 


were: Ng = Np = 


555, 


T-o = 0.2 


v = 


3.6, and p 


iiax 


= 25C 


A. 




































Faddeev 








Variational 




Adiabatic 


Potential 






equations 








methods 




approaches 


tin ax 


This work 


P 


M 


m @ 





1^0 




12| |3 


^ 




P 


I 


a 1 







0.084"' 0.0823 


0.082 0.092 








0.098 


HFDHE2 


2 


0.114"' 0.1124 


0.107 0.11 














4 


0.1167 






0.1171 


0.1173 













0.096"' 0.0942 


0.096 














HFD-B 


2 


0.131"' 0.1277 


0.130 
















4 


0.1325 






0.1330 




0.1193 0.133 0.131 


0.129 









0.0891 














0.106 


LM2M2 


2 

4 


0.1213 
0.1259 






0.1264 








0.1252 







0.0890 
















TTY 


2 

4 


0.1212 
0.1258 






0.1264 




0.126 







"'Results from ^ for a grid with Ng = Np = 275 and pn 



:60A. 



r(i) 



TABLE III. Excited state energy Ej. results for the helium trimer. The (absolute) values of E^ are given in mK. The grid 



(1) 



parameters 


used 


were: Ng - 


--Np- 


= 805, TO = 0.2, uo 


= 3.6, 


and pmax 


= 300 


A 




















Potential 




tmax 




This work 


I 


L9 


1 I 


Ti 


1 1 


21 


1 I 


16| 


117 


1 











1.5"' 1.46 


1 


46 


1.0' 


1 


1.517 




HFDHE2 




2 

4 




1.7"' 1.65 
1.67 


1 


.6 




1.665 















2.5"' 2.45 














HFD-B 




2 

4 




2.8"' 2.71 
2.74 








2.734 















2.02 










2.118 




LM2M2 




2 

4 




2.25 
2.28 








2.271 




2.269 











2.02 














TTY 




2 

4 




2.25 
2.28 








2.280 







"'Results from [Q for a grid with Ng = Np = 252 and pmax 



250 A. 
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TABLE IV. Trimcr excited-state energy E^ (niK) obtained with the TTY potential for various grids. 



^niax 


Ng^Np^ 252 
Pmax = 250'A 


Ne=Np = 502 

Pmax = 300 A 


A^e = iVp = 652 

Pmax = 300A 


Ng ^ Np = 805 

Pmax = 300 A 


Ng=Np^ 1005 

Pmax = 300 A 



2 


-2.108 
-2.348 


-2.039 
-2.273 


-2.029 
-2.258 


-2.024 
-2.253 


-2.021 
-2.248 



TABLE V. Estimations for *He atom-'*He dimer scattering lengths ^sc and inverse wave numbers >c ^ corresponding to the 
excited-state energy eI ' for the HFD-B, LM2M2 and TTY potentials. The accuracy for the scattering lengths is within ±5 A. 



The grid parameters used for the calculation 


of 4c are: Ng = Np ^ 502, tq = 0.18, f = 3.45 and p^ax 


= 460 A. 




Potential 


tmax -t-sc (A) 


X-' (A) 


Potential 


tmax 


4c (A) 


--' (A) 


HFD-B 


170"' 168 
2 145"' 138 
4 135 


109 
94 
93 


LM2M2/TTY 



2 
4 


168 
134 
131 


113 
98 
96 



""'Results from Q for a grid with Ng = Np = 320 and pmax = 400 A. 



TABLE VL Phase shift 5q '"' results (in degrees) for the HFD-B potential for various cm. energies E (in mK). The grid 
parameters used are: Ng = Np = 502, to — 0.18, ly — 3.45, and pmax =460 A. 



E 


"o 






E 








E 


"o 




"o 


-1.68541 


359.9 


359.9 


359.9 


-1.05 


299.1 


308.2 


309.2 


0.95 


262.4 


272.1 


273.7 


-1.68 


352.6 


353.9 


354.1 


-0.8 


290.8 


300.4 


301.5 


1.2 


260.0 


269.6 


270.7 


-1.65 


341.7 


345.0 


345.4 


-0.55 


284.4 


294.2 


295.4 


1.45 


257.8 


267.3 


268.4 


-1.60 


330.8 


337.7 


338.2 


-0.3 


279.3 


289.3 


290.4 


1.7 


255.9 


265.2 


266.3 


-1.55 


326.9 


332.8 


333.5 


-0.05 


275.1 


285.2 


286.3 


1.95 


254.1 


263.4 


264.5 


-1.50 


322.4 


329.0 


329.8 


0.2 


271.4 


281.3 


282.5 


2.2 


252.5 


261.7 


262.7 


-1.40 


315.4 


323.0 


323.9 


0.45 


268.1 


277.9 


279.0 


2.45 


251.0 


260.1 


261.1 


-1.30 


309.9 


318.1 


319.1 


0.7 


265.1 


274.8 


276.0 











TABLE 
Table Ivil 



VH. Phase shift 5, 



a„ 



k) 



results for the LM2M2 potential. The units and grid parameters used are the same as in 



E 




"o 


E 


"o 


gi:2) 


E 


^0 


4^' 


-1.30348 


359.8 


359.9 


-0.8 


304.6 


313.8 


0.95 


267.0 


276.2 


-1.3 


354.1 


355.3 


-0.55 


295.2 


304.8 


1.2 


264.1 


273.2 


-1.25 


337.9 


342.3 


-0.3 


287.9 


297.7 


1.45 


261.5 


270.6 


-1.20 


330.5 


336.3 


-0.05 


282.3 


292.2 


1.7 


259.2 


268.1 


-1.15 


325.2 


332.0 


0.2 


277.7 


287.4 


1.95 


257.1 


266.0 


-1.10 


321.1 


328.5 


0.45 


273.7 


283.2 


2.2 


255.3 


264.0 


-1.05 


317.6 


325.5 


0.7 


270.1 


279.5 


2.45 


253.6 


262.3 
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TABLE VIII. Phase shift 4'" 
0. 


"'"'' results for the TTY potential. The units and grid parameters used are 


the same as in Table 


E 








E 




^(2) 


5(4) 


E 




^(2) 




-1.30961 


359.7 


359.8 


359.8 


-0.8 


304.3 


313.5 


314.6 


0.95 


266.8 


276.1 


277.2 


-1.308 


355.9 


356.8 


356.9 


-0.55 


295.0 


304.6 


305.7 


1.2 


264.0 


273.1 


274.2 


-1.3 


350.2 


352.1 


352.4 


-0.3 


287.7 


297.5 


298.7 


1.45 


261.4 


270.5 


271.5 


-1.25 


336.8 


341.4 


341.9 


-0.05 


282.0 


292.0 


293.2 


1.7 


259.1 


268.1 


269.1 


-1.2 


329.7 


335.7 


336.4 


0.2 


277.5 


287.3 


288.4 


1.95 


257.0 


265.9 


266.9 


-1.10 


320.5 


328.1 


329.0 


0.45 


273.5 


283.1 


284.2 


2.2 


255.0 


263.9 


265.0 


-1.05 


317.1 


325.1 


326.1 


0.7 


270.0 


279.4 


280.5 


2.45 


253.5 


262.2 


263.2 



TABLE IX. The parameters for the *He— ''He Aziz and co-workers potentials used. 



Parameter 



HFDHE2 [^ 



HFD-B 



LM2M2 §1 ' 



£(K) 

rm (A) 

A 

a 

P 

Ce 

Cs 

Cio 

D 

Aa 

Cl 

C2 



10.8 
2.9673 

544850.4 

13.353384 



1.3732412 

0.4253785 

0.178100 

1.241314 



10.948 

2.963 

184431.01 

10.43329537 

-2.27965105 

1.36745214 

0.42123807 

0.17473318 

1.4826 



10.97 

2.9695 

189635.353 

10.70203539 

-1.90740649 

1.34687065 

0.41308398 

0.17060159 

1.4088 

0.0026 

1.003535949 

1.454790369 



TABLE X. The parameters for the "He-^He TTY potential used. 



A(K) 


315766.2067") 


Ce 


1.461 


P ((a.u.)-^) 


1.3443 


Cs 


14.11 


D 


7.449 


Cio 


183.5 


N 


12 







"'The value of A was obtained from the data presented in [BJ using, for converting the energy units, the factor 1 K= 3.1669 X 10 ^' : 



13 



^0, degrees 



360 



340 



320 



300 



280 



260 




-Blab, mK 
FIG. 1. S-wave helium atom - helium dimer scattering phase shifts (5o(-Eiab), i^iab ~ f (£- + kd|), for the HFD-B and TTY 



He- He potentials. The lower curve corresponds to the case where In 



while for the upper In 
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